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Abstract 

The importance of a Coulomb correction to the formalism proposed by 
Albergo et al. for determining the temperatures of nuclear systems at break- 
up and the densities of free nucleon gases is discussed. While the proposed 
correction has no effect on the temperatures extracted based on double isotope 
ratios, it becomes non-negligible when such temperatures or densities of free 
nucleon gases are extracted based on multiplicities of heavier fragments of 
different atomic numbers. 

PACS numbers: 25.70.Mn, 25.70.Pq, 25.70.-z, 24.10.Pa 

The formalism by Albergo, Costa, Costanzo, and Rubbino (ACCR) QTJ] offers simple and 
elegant prescriptions for the experimental determination of the temperature T, and of the 
free nucleon densities, p n p and p v p in a nuclear system at the instance of the break-up. 
This formalism presumes a specific scenario for the decay of excited systems, similar to 
the scenarios modeled by the Berlin Microcanonical Metropolis Monte Carlo || and the 
Copenhagen Statistical Multifragmentation || models. Therefore, it may be expected to be 
meaningful in circumstances (in an excitation energy domain), where the use of the above 
two more complete models can be justified. On the other hand, the ACCR formalism is 
incompatible with models such as the equilibrium-statistical model GEMINI |4]] and the 



Expanding Emitting Source Model ||, both of which refer to sequential decay scenarios of 
systems of uniform density. The purpose of the present paper is to point out that in certain 
circumstances, the approximation of the free nucleon g clS clS db collection of non-interacting 
nucleons, as assumed in the original ACCR approach, may not be sufficiently accurate, and 
that the inclusion of a proper Coulomb correction term is warranted. Regardless of the 
magnitude of the effects of such a correction, its inclusion is warranted already by didactical 
considerations. 

The ACCR formalism Q refers to a fragment production scenario in which an equili- 
brated freeze-out /break-up configuration emerges from initial compression and expansion 
stages. In this formalism, the average numbers of fragments or clusters of different mass and 
atomic numbers (A, Z) are determined by the requirement of chemical equilibrium between 
the fragments and free nucleons - neutrons and protons (quite obviously, the fragments are 
then also in a state of chemical equilibrium among themselves). In Ref. 1, the chemical equi- 
librium between fragments (A, Z) and gases of free neutrons and free protons is described 
by the equation 

P>A,z ~ B A ,z = Zji pF + (A - Z)p, nF , (1) 

where B A ,z is the fragment binding energy (taken with positive sign) and £ia,z, finF, and 
fl p F are the reduced chemical potentials of a fragment (A, Z), of a free neutron, and a free 
proton, respectively. Here, the qualifier "reduced" is used with respect to the term chemical 
potential and a tilda is used in the respective symbolic representation, to distinguish the 
quantities involved in Eq. [1] from the true thermodynamical chemical potentials Ha,Zi as 
defined via the equation 

_ dF(V,T) 

»a,z - [-qn^\v,t- (2) 

In Eq. |2], F(V, T) is the free energy of the system at volume V and temperature T. The 
quantity N A .z is the average number of fragments (A, Z) , and the partial derivative is taken 
at constant temperature and volume. Note that, unlike their "reduced" counterparts p>A,z, 
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the true chemical potentials fi A ,z include not only the binding energy term Ba,z, but also the 
energy of interaction of the fragments (A, Z) with the relevant mean Coulomb field. When 
true chemical potentials are considered, the chemical equilibrium is expressed through the 
following equation: 

fi At z = Zfj, pF + (A - Z)^ nF - v£rt{Z), (3) 

where V^ nU (Z) is the average energy of mutual Coulomb interaction of Z free protons. The 
latter Coulomb interaction energy must be subtracted on the right-hand side of Eq. ||| to 
compensate for the fact, that in a mean- field type of "bookkeeping" of the Coulomb energy 
that is included in the term Zfi pF , the mutual Coulomb interaction energy is double-counted. 

For some purposes, one may consider the convenient notion of a reduced chemical po- 
tential jl(A, Z), which excludes the fragment binding energies and the interaction with the 
mean Coulomb field. In Maxwell-Boltzmann statistics, the latter quantity has a simple re- 
lationship |]J to the density p A ,z of fragments (A, Z) representing the number of fragments 
per unit break-up volume: 



= A^m^ zlT (4) 

An 



Eqs. |] and Eq. p] are the two fundamental equations used in Ref. [1] to establish the 
relationship between the yields of various fragments and the characteristics of a chemically 
equilibrated nuclear system at break-up. In Eq. f|, At = h/y/2irm T is the nucleon ther- 
mal wave-length (m is the mass of a nucleon), and uj a ,z(T) is the temperature-dependent 
internal partition function of fragment (A, Z): 

^,z(T) = S fe (2^ + l)e-^' Z / r , (5) 

where the summation extends over all bound states of the fragment (A, Z) with spins s k ' 
and excitation energies E^' z . 

An inspection and comparison of Eq. [TJ, to the more fundamental Eq. |^ reveals a lack of 
symmetry of the former equation. While the fragment (A, Z) side of the balance includes 



the mutual interaction energy of the constituent nucleons - the binding energy Ba,z, no 
equivalent term is present for the free nucleons on the r.h.s. of Eq. [1]. Yet, the Z free 
protons do interact among themselves via long-range Coulomb interactions. Therefore, 
a more complete equation for the chemical equilibrium based on Eq. |3| must include the 
respective Coulomb interaction term. Note that, on the fragment side of the balance, the 
mutual Coulomb interaction energy of Z protons is included in the binding energy term 
Ba,z- One may note also that neither side of Eq. |l| considers explicitly (or implicitly) the 
Coulomb interaction energy of the Z protons with the remaining (Z syatem — Z) "spectator" 
protons of the system. Such an omission, however, may be well justified, as these two 
Coulomb interaction energies are to a good approximation equal to each other and, hence, 
cancel mutually. 

A more complete, Coulomb-corrected equation for the chemical equilibrium of fragments 
(A, Z) and free nucleons, in terms of reduced chemical potentials pLA,z has the symmetrical 
form: 

P-a,z ~ B A ,z = Zfi pF + (A - Z)fl nF + V^JZ). (6) 

Here, V^{Z) represents the average potential energy of the mutual Coulomb interaction 
of Z free protons, an equivalent of the term -Ba,z- It is worth noting that, here (unlike in 
Eq. HD the Coulomb interaction term enters with positive sign, as no Coulomb interaction is 
included in the reduced free-proton chemical potential p, p p. 

While it is clear from simple estimates that the Coulomb correction term, VjJ^^Z), is of 
non-negligible magnitude when compared to typical temperatures of the system, it is not 
obvious how to actually evaluate it. A conservative estimate for the value of V^(Z) may 
be obtained by assuming that this term is equal to the Coulomb interaction energy of Z 
protons uniformly distributed over a spherical volume Vf ree = Z/p pF , i.e., distributed with a 
density equal to that of the gas of free protons p pF at break-up. In this case, the correction 
term is independent of the fragment mass number A: 
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V cLm = \ - f z Z l /3 « l.39Z^p^(MeV), (7) 

where p p i? is expressed in units of fm~ 3 . 

The presence of the Coulomb correction term V^(Z) in Eq. |] modifies the basic equation 
5 of Ref. [1] for the average number of fragments (A, Z) per unit break-up volume, pa,z- It 
can now be written more accurately as 

Pa,z = ^^^^ ^^^^^^/^ (8) 

replacing Eq. [|. In Eq. [§], p n p and p p p are the densities (i.e., numbers per unit break-up 
volume) of free neutrons and free protons, respectively. A similar result was obtained earlier 
H based on a more rigorous macrocanonical description of a decaying nuclear system in a 
freeze-out configuration. It is worth noting that, in more complete theoretical descriptions of 
equilibrated freeze-out configurations, offered by the Berlin and Copenhagen j|] models, 
the effects of the Coulomb interaction of the free protons are accounted for in a rigorous 
fashion, but remain largely transparent to the model users. 

In practical applications of Eq. [|, ratios of properly selected densities, Pa 1 ,z 1 / Pa 2 ,z 2 , 
are taken and identified with the ratios of the respective experimental yields of fragments 
Ya 1 ,z 1 /Ya 2 ,z 2 - Such ratios are free of some model parameters (e.g., of the densities of free 
neutron and proton gases, in the case of double isotope ratios), providing often a simple link 
between observable yields and selected characteristics of the break-up state. 

It is clear from Eq. ||that the introduction of the Coulomb correction term V^i(Z) is of 
no consequence when ratios of yields are taken for fragments with identical atomic numbers 
Z, i.e., ratios of experimental fragment yields of the type Ya+i,z/Ya,z- In such cases, the 
corresponding Coulomb correction terms for the two isotopes involved cancel each other. 
As a result, this correction has no effect on the outcome of an experimental evaluation 
of break-up temperatures based on double-isotope ratios - the most common use of the 
ACCR UH approach. A similar cancellation does not, however, occur in cases when, e.g., an 
experimental "thermometer" is constructed from isotone ratios, Ya+i,z+i/Ya,z, or in cases 



when relative densities of free neutron and proton gases, p n F/p P F, are determined based on 
an isobaric ratio Y^,z /Ya,z+i- To assess the significance of the proposed Coulomb correction, 
several examples are considered below. In these examples, it is assumed that T=3.3 MeV 
(as found for the system S+Ag at E/A=22 MeV) and p pF = 5/(4/3tt8.0 * 141) = 0.0011 
fm -3 (which corresponds to 5 protons in a break-up volume of radius Rbreak-u P = 2.0* 141 1 / 3 
fm, as for the system 32 S+ 109 Ag). 

First, consider an evaluation of the ratio of densities of free neutron and free proton gases 
from the observed isobaric ratios Y^z/Y^z+i- Based on Eq. ||, one has 

P nF _ p ztWp /q\ 

— flraw-r Couh \ J I 

PpF 

where R raw is the value of this ratio deduced in the absence of the Coulomb correction (i.e., 
given by the original ACCR ]1] formalism) and F^ ul is a correction factor resulting from 
the Coulomb term proposed in the present paper: 

F n J: ul = ett^ 1 )-^ W (10) 

Using Eqs. |Tt] and [7| and the values of T=3.3 MeV and p p p = 0.0011 fm -3 , one obtains 
Fcoui = 11, for the case of the isobaric ratio Y^i/Y 32 (tritium - helium-3), FqJ^ = 1.28, for 
the case of the isobaric ratio V 13 6 /Y" 13 7 , and FqJ^ = 1.59 when the isobaric ratio V 3416 / Kj4 il7 
is utilized. This example demonstrates that even in the favorable case of light isobars 3 if and 
3 He, the correction factor is large enough to mandate an inclusion of the proposed Coulomb 
correction term in the equation for the chemical equilibrium. Certainly, this correction factor 
is quite sizeable when yields of heavier isobars are utilized for the evaluation of the relative 
densities of free neutron and free proton gases. 

It is worth noting that, according to Eq. [10], the Coulomb correction factor, F^ ul) for 
the relative densities of gases of free neutrons and protons is always greater than unity, since 
the exponent [V^(Z + 1) — V^(Z)] is positive. This fact reflects the role of the Coulomb 
energy in a neutron-enrichment of the free nucleon gas, a rather trivial effect that should not 



be confused with an isospin fractionation |8]-|r0"[ driven by an isospin-dependent equation of 
state of nuclear matter. 



As a second example, consider the evaluation of the break-up temperature T based on a 
double isotone ratio: 

-^isotone \r Ixr ' \ ) 

Y M,Z 2 I Y A 2 +1,Z 2 +1 

In such a case, the breakup temperature T is ultimately evaluated from the experimen- 
tally determined value of the ratio 



AB + AV c p f ul AB 



ln(aR isotone ), (12) 



T T 

- 1 - 1 raw 

where AB = B Al+hZl+1 -B AuZl +B A2jZ2 -B A2+h z 2+ i and AV p Q F ul = Kw(^i + l)-V^(Zi) + 
Vcoui(Z2) — V£oui(Z2 + 1)- The parameter a in Eq. [TJ] accounts for intrinsic partition functions 
uj (see Eq. 0) of the isotones involved, and T raw is the break-up temperature obtained using 
the original ACCR [|l|] approach. 

Eq. [12] allows one to express the relevant Coulomb correction factor Fq ou1 as 

T AV pF 
Fl ul = ^- = 1 + ^f. (13) 

l raw L±D 

Using Eq. the value of p v p = .0011 fm -3 , and values of binding energies from 
the mass tables, one obtains Fq ou1 = 0.96 in the case when the experimental double 
isotone ratio (Yi 3 j/Yi2fi)/(Y^2/Y 3t i) is employed, Fq ou1 = 1.19 in the case of the iso- 
tone ratio (Yi 3 j/Yi2fi)/(Y Wt5 /Y 9 ^), and F£ oul = 0.31 in the case of the isotone ratio 
(^14,8/^13,7)/ (^3, 2/^2,1)- Again, the estimated magnitude of the Coulomb correction fac- 
tor well warrants an inclusion of the Coulomb correction term in the ACCR formalism. 

In summary, the importance of a Coulomb correction term to the equation for the chem- 
ical equilibrium between fragments and the gas of free nucleons has been demonstrated. 
The correction term restores the symmetry of the equation defining the equilibrium, when 
the mutual interaction energies of nucleons in both, bound and free states, are consistently 
accounted for as done with the more rigorous Eq. [5| While the proposed correction term has 
no effect on the determination of break-up temperatures based on double isotope ratios and 
may be small in some cases, its effects on the determination of break-up temperature from 



isotone ratios and on the determination of the relative densities of free proton and neutron 
gases may be quite substantial in some other cases. This correction is certainly important in 
systematic studies of various experimental "thermometers" that, by design, include a large 
variety of isotonic ratios. 
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